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1. INTRODUCTION 

Surface reflectance retrieval from imaging spectrometer data as acquired with the 
Airborne Visible/In flared Imaging Spectrometer (AVIRIS) has become important for 
quantitative analysis. In order to calculate surface reflectance from remotely measured radiance, 
radiative transfer codes such as 5S (Tanr6 et al., 1990) and MODTRAN2 (Berk et al., 1989) 
play an increasing role for removal of scattering and absorption effects of the atmosphere. 
Accurate knowledge of the exo-atmospheric solar irradiance (EJ spectrum at the spectral 
resolution of the sensor is important for this purpose. The present study investigates the 
impact of differences in the solar irradiance function, as implemented in a modified version of 
5S (M5S) (Teillet and Santer, 1991), 6S (Vermote et al., 1994), and MODTRAN2 (Berk et 
al., 1989), and as proposed by Green and Gao (1993), on the surface reflectance retrieved from 
AVIRIS data. Reflectance measured in situ is used as a basis of comparison. 

2 . DATA USED 

An asphalt site within an AVIRIS scene acquired over the Greater Victoria Watershed 
near Victoria, British Columbia on August 29, 1993 was selected for this study. The 
extracted AVIRIS radiance spectrum was averaged over the target of 2 x 3 pixels, which 
excludes edge pixels. Corresponding ground-based reflectance data were acquired with a GER 
MARK V spec troradiometer and a Spectralon panel. These data are accurate to within ±5% 
reflectance. Solar attenuation measurements were collected on the ground at the time of the 
AVIRIS overflight with a SONOTEK sunphotometer for the calculation of atmospheric 
optical depths. In addition, estimated water vapour concentrations as reported by Staenz et al. 
(1994) for this particular AVIRIS dataset were used for atmospheric modelling. 

Several E 0 data sets were included in this study (Table 1). These E o spectra are based 
on Iqbal (1983) for the MSS code, Neckel and Labs (1984) for 6S, and a combination of 
Neckel and Labs, Wehrli (1985), and Thekeakera (1974) for MODTRAN2. The proposed 
update of the E 0 function in MODTRAN2 suggested by Green and Gao (1993) was also 
included. This E 0 spectrum is based on Neckel and Labs (1984) and data provided by the 
ATMOS sensor onboard the Space Shuttle. 

3. ANALYSIS AND RESULTS 

Three radiative transfer codes, MSS, 6S and MODTRAN2, were used in combination 
with each of the four aforementiond Eq functions to generate a total of 12 computed surface 
reflectances for the asphalt site. The procedure was carried out in three steps. First, the 
radiative transfer codes were used to generate tables of simulated AVIRIS radiance spectra. 

Then, the tables were scaled by the spectral ratio of the Eo function under study to that which 
is implemented in the radiative transfer code (both convolved to the AVIRIS handset). 

Finally, multiplicative and additive coefficients for atmospheric correction were computed 
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from the tables and applied to the airborne radiances (Williams et al., 1992). 

The Eo functions are plotted for the AVIRIS bands in Figure 1. Variations of up to 
20% occur between the different data sets, although the differences are generally within ±5%. 
The main differences are found in the 2100 - 2500 nm region and other notable discrepancies 
are apparent around 940 nm, 1300 nm, 1650 nm, and 1900 nm. The 6S Eo function follows 
the Green and Gao function very closely with the exception of the 1300 nm region. A similar 
trend can generally be seen for the M5S and the MODTRAN2 E 0 spectra above 800 nm, 
although deviations from this tendency occur for some of the AVIRIS bands above 1800 nm. 

These dev iati ons between the E 0 functions translated to similar relative differences in 
the retrieved surface reflectances for each of the radiative transfer codes used. Figure 2 shows 
the reflectances obtained with the MODTRAN2 code. Figure 3 is a blow-up of the 2050 - 
2450 nm region of Figure 2. The reflectance spectra generated with the 6S and Green and Gao 
Eo functions are slightly smoother in this wavelength region than those produced with the 
original MODTRAN2 or M5S Eo datasets. However, none of the retrieved reflectances exactly 
matches the ground-based spectrum. This is also true for the 940 nm, 1300nm, and 1600 nm 
regions (Figure 2), although a closer match with ground-based reflectances is apparent in the 
last two regions for the reflectances generated with the 6S/Green and Gao Eo functions. 

Changes in the overall shape of the surface reflectance spectra due to the different Eo 
functions (Figure 3) are less than those generated by the various radiative transfer codes using 
the same E 0 dataset (Figure 4). Differences arising from the use of different radiative transfer 
codes occur especially in the wavelength regions affected by atmospheric gases. For example. 
Figure 4 shows major reflectance differences above 2250 nm between MSS and 
6S/MODTRAN2 generated spectra, even though the Eo function used (Green and Gao) was 
the same in each case. This is mainly due to the different methods used for calculation of the 
atmospheric transmission in these codes (Staenz et al., 1994). 

4. CONCLUSIONS 

Differences in the E 0 spectra as implemented in M5S, 6S and MODTRAN2 radiative 
transfer codes and as proposed by Green and Gao can translate to relative differences in retrieved 
surface reflectance of up to 20%, although they are generally within ±5% at most 
wavelengths. The largest deviations are found above 21 00 nm. In addition, notable 
differences occur in the 940 nm, 1300 nm, and 1650 nm regions. The overall shape of the 
retrieved reflectance spectrum is generally less affected by differences between Eo functions 
than by differences between the atmospheric radiative tansfer codes. 
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Table 1: Exo-atmospheric solar irradiance functions (E c ) and their spectral resolution for the AVIRIS 
wavelength coverage (400-2500 nm) with respect to different radiative transfer codes. 


Source 

Radiative Transfer Code 
(spectral resolution) 

Spectral Resolution of E 0 Function H 

original 

as implemented in radiative 
transfer code 

Iqbal 

MSS (5 nm, 20 cm*0 

5-100 nm 1 

5 nm 

1 Neckel and Labs 

6S (2.5 nm) 

1-5 nm 2 

5 nm 

Neckel and Labs, 

Wehrli, 

Thekeakera 

MODTRAN2 
(variable, finest = 1 cm 1 ) 

1-5 nm 2 
1-10 nm 3 
5- 100 nm 1 

20 cm* 1 

Green and Gao 

proposed for MODTRAN 

2.5 nm 


1 5nm (400-6 lOnm), 

10 nm (610-1000 nm), 50 nm 

(1000-2000 nm), 100 nn 

n (2000-2500 nm) 


2 I nm (400-630nm), 2 nm (630*869 nm), 5 nm (869-1248 nm) 

3 1 nm (400-630nm), 2 nm (630-999 nm), 5 nm (999-2002 nm), 10 nm (2002-2500 nm) 
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Percent differences between the exo-astmospheric solar irradiance (Eo) spectra with respect 
to the Green and Gao Eo spectrum. The irradiance data sets were convolved to the 
AVIRIS bands. 
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Surface reflectances of asphalt retrieved from AVIRIS data with MODTRAN2 using 
different solar irradiance functions, and the corresponding ground-based reflectance (GER). 
The atmospheric conditions used in the retrieval were: water vapour = 1.5 g/cm 2 ; aerosol 
optical depth at 550 nm = 0.121. 



Surface reflectances of asphalt in the 2050-2450 nm region (Blow-up of Figure 2). 



Surface reflectances of asphalt retrieved from the AVIRIS data with the different radiative 
transfer codes using the solar irradiance function proposed by Green and Gao as well as the 
corresponding ground-based reflectance (GER). The models were run with the same 
atmospheric conditions as those noted in Figure 2. 
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